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Abstract. We consider a steady state Vq of the Euler equation in a fixed bounded domain 
in R". Suppose the linearized Euler equation has an exponential dichotomy of unstable 
and center-stable subspaces. By rewriting the Euler equation as an ODE on an infinite 
dimensional manifold of volume preserving maps in W k ' q , (k > 1 + the unstable (and 
stable) manifolds of vq are constructed under certain spectral gap condition which is verified 
, for both 2D and 3D examples. In particular, when the unstable subspace is finite dimen- 

£^ ■ sional, this implies the nonlinear instability of vo in the sense that arbitrarily small W k,q 

perturbations can lead to I? growth of the nonlinear solutions. 
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1. Introduction 



We consider the incompressible Euler equation on a smooth bounded domain Q GG 
n > 2, under the slip (or periodic in certain directions) boundary condition 

v t + (v ■ V)i> = — Vp and V • v = x G Q 

v ■ N = 0, xedVl 



(E) 



where v — (v 1 , . . . , v n ) T is the velocity field and N is the unit outward normal vector of Q. 



(■Sj | We take 

in 



(1.1) W^ ler = {v GW k ' q {Q,R n ) | V-v = 0intt, v-N = 0ondQ}, q > 1, k > 1 + ^ 



as the phase space. It is well known that (E) is well posed in these spaces, globally if n = 2 
and locally if n > 3. As shown below, the pressure p can be written in terms of v through a 
quadratic mapping. 

Let vo be a steady solution of (E). Linearize the equation at (E) and we obtain 
(1.2) v t = -{v ■ V)u - (v ■ V)v - Vp = Lv 

where the operator L can be defined as acting only on v since the linearized pressure p 
can be determined by v linearly, though non-locally. To study the dynamics near vq, the 
first step is to understand linear instability, that is, the spectrum of the operator L. The 
problem of linear instability of inviscid flows has a long history dated back to Rayleigh and 
Kelvin in 19th century. But even until now, very few sufficient conditions for the existence of 
unstable eigenvalues are known and most of the investigations had been restricted to shear 
flows and rotating flows. See [DR81J,|DH66] and the references therein. Some recent results 
on instability conditions can be found in [L03] |L05] for shear flows and rotating flows, and 
in |L04aj for general 2D flows. Besides the discrete unstable spectrum, the linearized Euler 
operator may also have non-empty unstable essential spectrum due to nontrivial Lyapunov 
exponents of the steady flow v ( [FV91j [LM91j [SL09] [LS03] jV96] b Indeed, growth of 
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linearized solutions can be seen in H s (s > 1) norm near any nontrivial steady flows, due to 
the stretching of the steady fluid trajectory. One also notes that the choice of the Sobolev 
space (norm) actually affects the essential spectrum which corresponds to small spatial scales, 
but not discrete spectrum corresponding to large scales. 

Consider a linearly unstable steady flow vq, that is, the linearized Euler operator L has an 
unstable discrete eigenvalue. To discuss the nonlinear instability, it is important to specify 
the norms. On the one hand, certain regularity is necessary in the local well-posedness 
of classical solutions. On the other hand, as mentioned in the above, the choice of the 
norm already affects the essential spectrum at the linear level. Moreover, even near steady 
states without unstable eigenvalues, solutions are expected to grow in H s norm with s > 1. 
Therefore the growth in the energy norm L 2 of nonlinear solutions is more a nonlinear 
reflection of the linear instability from the discrete spectrum, which also corresponds to the 
instability in the large scale spatial scale (see [L04bt Section 6.2] for detailed discussions). 
Naturally, the ideal nonlinear instability result would be to obtain order 0(1) growth in L 2 
distance (weaker energy norm) from the steady state vq of solutions starting with arbitrary 
small initial perturbation from vq in H s norm (stronger norm), where s > is determined by 
the regularity of unstable eigenfunctions. Such nonlinear instability result (i.e. H s — > L 2 ) is 
not only mathematically stronger, but also physically interesting due to the discussions in 
the above. 

The proof of the nonlinear instability based on unstable eigenvalues is nontrivial for several 
reasons. The main difficulty is that the nonlinear term v ■ Vv contains a loss of derivative. 
Moreover, the norm- dependent unstable essential spectrum corresponds to growth in small 
spatial scales. It may interact with discrete unstable modes and then cause complications in 
proving nonlinear instability. In the last decade, there appeared several proofs of nonlinear 
instability for Euler equations f |BGS02j [GreOOj [VF03] [FSV97] [L04bj). In [FSV97j . non- 
linear instability in H s (s > 1 + |) norm was proven for n— dimensional Euler equations, 
under a spectral gap condition which was verified for 2D shear flows. All other papers proved 
nonlinear instability only for the 2D case, in the more desirable H 1 or L 2 norms and under 
different spectral assumptions. Particularly, in |L04bj nonlinear instability in L 2 norm was 
proved for general linearly unstable flows of 2D Euler, without additional assumption on the 
growth rate which was made in QBGS02] jGreOO] |VF03] ). 

When there exist a collection a u of unstable eigenvalues of the linearized Euler operator 
at the steady flow vq with strictly larger real parts than the rest of the spectrum, it is 
very natural from a dynamical system point of view to ask if a locally invariant unstable 
manifold tangent to the eigenspace of a u exists. The answer to this question would provide 
a better picture of the nonlinear instability, including the dimensions and the directions 
of the unstable solutions with certain minimal growth rate. More importantly, such locally 
invariant manifolds provide more precise characterization of the local dynamical pictures near 
an unstable equilibria, and are also basic tools for constructing globally invariant structures 
such as heteroclinic and homoclinic orbits. These dynamical structures are important in 
understanding the turbulent fluid behaviors. The major obstacle to the construction of 
local invariant manifolds is again the loss of derivative due to the derivative nonlinearity 
v-Vv. For dissipative models such as reaction-diffusion equations ( |H81j ) and Navier-Stokes 
equations f |Yu89] |Li05j ). it is rather standard to construct invariant manifolds since the 
dissipation terms provide strong smoothing effect to overcome the loss of derivatives in the 
nonlinear terms. However, for non-dissipative continuum models including Euler equations, 
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the linearized operators have no smoothing effect to help overcome the loss of derivative 
in the nonlinear terms. So it has been largely open to construct invariant manifolds for 
conservative continuum models such as Euler equations. In the proof of nonlinear instability 
for 2D Euler, one takes initial data perturbed along the direction of unstable eigenfunctions 
and then uses special properties of nonlinear solutions of 2D Euler to overcome the loss of 
derivatives, such as the bootstrap arguments in ([BGS02] |L04b] |VF03J) and the nonlinear 
energy estimates in ( |Gre00j ). However, these techniques can not be used for constructing 
invariant manifolds, since we do not know beforehand the initial conditions for solutions on 
unstable or stable manifolds. 

In this paper, we obtain the first result on stable and unstable manifolds of Euler equations 
in any dimensions. To state the precise results, we formulate the following assumptions: 

(Al) v G W^, where r > 4. 

(A2) 3 X u > A cs > 0, and closed subspaces X u and X cs of W^ ler such that they satisfy 
L(X U ^ CS fl domain(L)) C X UtCS respectively and W E '£ ler = X u © X cs . Moreover, let 
L u = L\x u and L cs = L\x cg , then for some M > 0, they satisfy 

\e tLcs \ < Me Xcs \ V t > and \e tLu \ < Me Xu \ V t < 0. 

(A3) The largest Lyapunov exponent fio (in both forward and back time) of the linearized 
equations 

y t = Dv (x(t))y y t = -(Dv (x(t))} y, 

along any integral curve x(t) of x t = v (x), satisfies 

A u — A cs > K {1 

where Kq is a constant depending only on r and k. 
Now we give our main theorem 

Theorem 1.1. Under above assumptions (A1)-(A3), there exists a unique C r ~ 3,1 local un- 
stable manifold W u of v in W E '£ ler which satisfies 

(1) It is tangent to X u at vq. 

(2) It can be written as the graph of a C r ~ 3,1 mapping from a neighborhood of Vq in X u 
to X cs . 

(3) It is locally invariant under the flow of the Euler equation (E), i.e. solutions starting 
on W u can only leave W u through its boundary. 

(4) Solutions starting on W u converges to Vq at the rate e xt as t — >• — oo for any A < 

The same results hold for local stable manifold of vq as the Euler equation (E) is time- 
reversible. 

Remark. In Sections 3 and 4, the assumptions (A2)-(A3) are verified for linearly unstable 
2D shear flows and rotating flows, as well as 3D shear flows. For these flows, /i = and 
thus (A3) is automatically satisfied if an unstable eigenvalue exists. 

Remark. Suppose v E W hl+r,qi fl W k2+r ' q ' 2 and the invariant decompositions W ki ' Qi = X l u © 
X l cs along with the same exponents \ UfiS satisfy (A2) - (A3) fori = 1,2. One may construct 
the local unstable manifolds Wf C W ki,Qi , i = 1,2 from Theorem \l.l\ Assume X^ = X^, 
we claim = W% on an open neighborhood of vq. In fact, let k = max-f/ci,/^} and 
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q = max{<7i, 52} and W u be the unstable manifold of Vq in W k ' q C W ki,qi , i — 1,2. Clearly 
W u C , i — 1,2, with the same tangent spaces and thus the above claim follows. 

When X u is finite dimensional, which in particular is always true in 2D under assumption 
the \ u > A cs > [k — 1) /io as proved in Section[3], then the W k,q topology and L 2 topology are 
equivalent on W u . Then an immediate consequence of the above theorem is the nonlinear 
instability in L 2 norm with initial data slightly perturbed from vq in W h,q norm. 

Corollary 1. Suppose (Al) -(A3) are satisfied andX u is finite dimensional, then there exists 
5 > such that there exist a solution v(t) such that \v(0) — vq\ L 2 > 5 and \v(t) — vo\ W k, q — > 
as t — > —00 exponentially. 

This is a stronger statement than the usual exponential nonlinear instability which by 
definition means that there exists 5 > such that for any e > 0, there exists a solution 
v(t) satisfying both \v (0) — i> | < e and sup 0<t<O( -_ log(: ) \v(t) — v \ > 5. One notes that 
the perturbed solution v(t) is allowed to depend on e and no condition is imposed on the 
asymptotic behavior of v(t) as t — > —00. In the contrast, any solution v(t) in Corollary [T] 
satisfies, in addition to the requirements in the nonlinear instability definition for all e > 0, 
that it starts at Vq when t = —00 and get out of the (^-neighborhood of Vq in L 2 norm. 

In the previous works (see references above) on nonlinear instability of the Euler equation, 
growing solutions have usually been found in the most unstable direction of the linearized 
equation (11.21) with roughly the maximal exponential growth rate. The above unstable 
manifold theorem actually provides solutions growing in other relatively weaker unstable di- 
rections. Though not necessary, it is easier to see this when X u is finite dimensional. In fact, 
on the finite dimensional locally invariant manifold W u , the Euler equation (E) becomes a 
smooth ODE and Vq is a hyperbolic unstable node. When L\x u has eigenvalues with different 
real parts, one may split X u into strongly unstable subspace X uu and weaker unstable sub- 
space X wu such that X u = X uu ® X wu . The standard invariant manifold theory implies there 
exist the locally invariant weakly unstable manifold W wu tangent to X wu and the locally 
invariant strongly unstable fibers W£ u with base point v G W wu and extend in the direction 
of X uu . Those solutions in W wu grow in the directions of X wu at a slower exponential rate. 
Moreover, the Hartman-Grobman theorem implies that a Holder homeomorphism on X u 
may transform the Euler equation restricted on W u into a linear ODE system. 

In Section 4, we construct linearly unstable 3D steady flow satisfying the assumptions in 
Theorem 11.11 By Corollary [U this implies nonlinear exponential instability in L 2 norm. To 
our knowledge, this is the first proof of nonlinear instability of 3D Euler equation. We note 
that the methods for proving nonlinear instability of 2D Euler cannot be applied to prove 
nonlinear instability for 3D Euler. For example, the bootstrap arguments in ([BGS02J [L04b] 
[VF03] ) strongly use the fact that vorticity is non-streching in 2D and therefore do not work 
in 3D due to the vorticity stretching effect. 

Below, we sketch the main ideas in the proof of Theorem 11.11 The main difficulty in 
constructing the unstable manifolds for the Euler equation lies in the fact that a derivative 
loss occurs in the nonlinear terms while the linearized flow does not have the smoothing 
property. We will prove Theorem 11.11 mainly by considering the Euler equation (E) in 
the Lagrangian coordinates. In a seminal paper [Ar66j, V. Arnold pointed out that the 
incompressible Euler equation can be viewed as the geodesic equation on the group of volume 
preserving diffeomorphisms. This point of view has been adopted and developed by several 



UNSTABLE MANIFOLDS 



5 



authors in their work on the Euler equations, such as |EM70j ISh85t IBr99l ISZ08al ISZ08bj to 
mention a few. 

On the one hand, the main advantage of this approach is that (E) on a fixed domain as 
in this paper becomes a smooth infinite dimensional ODE |EM70] on the tangent bundle of 
the Lie group Q of volume preserving diffeomorphisms of Q and thus the difficulty arising 
from the loss of regularity disappears. A side remark is that this is in contrast with the 
Euler equation with free boundaries [SZ08a, SZ08b] where the Riemannian curvature of the 
infinite dimensional manifolds of volume preserving manifolds are unbounded operators and 
the Euler equation can not be considered as infinite dimensional ODEs. As a clarification, by 
saying that the Euler equation on fixed domains defines an ODE, we mean that it corresponds 
to a vector field which is everywhere defined and smooth on an infinite dimensional manifold. 
In this sense, evolutionary PDEs involving unbounded operators such as heat equation or 
wave equation do not define infinite dimensional ODEs. 

On the other hand, in the Lagrangian coordinates, the steady state Vo generates a special 
geodesic Uo(t) which coincides with the integral curve starting at the identity map of the right 
invariant vector field on Q generated by vq. By carefully using local coordinates along this 
orbit generated by the group symmetry, we further transform the localized Euler equation 
into a weakly nonlinear non-autonomous ODE with linear exponential dichotomy. One 
notes that the multiplication on this Lie group Q - the composition between W k,q volume 
preserving maps - is continuous, but not smooth, see Proposition 12.21 The assumption 
v ^ W k+r,q ensures that the localization we choose possesses certain smoothness. Moreover, 
compared to the usual exponential dichotomy enjoyed by many ODEs and even PDEs (see 
for example [CL88J). the exponential dichotomy here has the defect that the angle between 
the associated invariant subspaces may not have a uniform positive lower bound in time due 
to the possible growth of the linearized ODE flow of the vector field vq. Moreover, this type 
of unwanted growth in t may also appear in the norms of the nonlinearities. Assumption 
(A3) is used to overcome this non-uniformity in t of the dichotomy in the construction of the 
unstable manifolds via the method based on the Lyapunov-Perron integral equations. Here 
we have to take the advantage of the fact that solutions on the unstable manifolds decay to 
Vq exponentially as t — > — oo. Since solutions on the center manifolds do not have similar 
decay properties, we still can not construct center manifolds of stead states. 



Remark. The assumption (A2) is the standard linear exponential dichotomy for constructing 
invariant manifolds. The extra gap assumption (A3) might be technical, but appears rather 
natural in our approach. In the proof, we change back and forth between Lagrangian and 
Eulerian coordinates, and each transformation induces a factor e mt in the estimates. The 
extra gap \ u — A cs > -Ko/^o guarantees that after all these transformations, the exponential 
dichotomy of unstable and central-stable parts still persists. 



The method introduced in this paper provides a general approach to construct unstable 
manifolds for many other continuum models in fluid and plasmas. In these models, the loss 
of derivative is also due to nonlinear terms from the material derivative. By working on 
Lagrangian coordinates, we can again overcome such a loss of derivative and the existence of 
unstable manifolds is conceivable with sufficient spectral gap. We are using this approach to 
construct unstable manifolds for density-dependent Euler equations and the Vlasov-Poisson 
system for collisionless plasmas. 
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2. Proof of Theorem 11.11 

Lagrangian coordinates is a standard tool in studying the Euler equation. In Subsection 
12.11 and 12. 2\ we will present the manifold structure of the set Q of the Lagrangian maps and 
the ODE nature of (E) on TQ. These general results have actually been proved even for 
Euler equations defined on Riemannian manifolds in [EM70J in a rather geometric language. 
However, we need to establish a more concrete framework along with more detailed estimates 
to be used in the construction of local invariant manifolds in Subsections 12.31 - 12.61 which is 
done in a more directly equation based manner in Subsections 12. II and 12. 21 In Subsection I2.3[ 
we rewrite in the Lagrangian coordinates the localized Euler equation in a neighborhood of 
the solution curve generated by Vq and in Subsection I2.4| the linear exponential dichotomy 
is given. The unstable integral manifold corresponding the linear exponential dichotomy 
is constructed in Subsection 12.51 and then finally the unstable manifold in the Eulerian 
coordinates is obtained in Subsection 12. 61 

Throughout this section, we will use K > as a generic constant depending only on r 
and k and C > only on n, r, k, q, vq. Both K and C may change from line to line. We will 
use D or V to denote the differentiation with respect to physical variables in Q and T> the 
Frechet differentiations in function spaces. 

2.1. Lagrangian coordinates and the Lie group of volume preserving maps. Let 

u(t, -) : Q — > fl be the Lagrangian coordinate map defined by 

(2.1) u (0,y)=y and u t (t, y) = v(t, u(t, y)). 

In particular, let u (t, y) be the Lagrangian map of the steady vector field v (x). Throughout 
the paper, we fix a constant /i > fio > such that 

(2.2) X u - X cs > K tM 

where fio is the Lyapunov exponent of vq. From the definition of the Lyapunov exponent, 
we have 

(2.3) \u (t,-)\ c i + \(u (t,-)y 1 \ c , <Ce lflltl , tER,0<l<k + r 

for some C > independent of t. This possible exponential growth of the norm of uq 
makes the problem much more subtle than the unusual constructions of the local invariant 
manifolds in differential equations. 

Since the flow is incompressible and k > 1 + -, we have for any t G R, 

(2.4) u(t, •) G G = {</> G W k ' q (tt,W n ) | (f) is a diffeomorphism, det = l,(j)(dQ) = dtt}. 

Clearly the composition makes Q a group. We will show that Q is an infinite dimensional 
submanifold of W/ fc ' 9 (fi,M n ). This will be our configuration space when the Euler equation 
is written in the Lagrangian coordinates. We will work with local coordinates on Q. 
Formally, the tangent space of Q is given by 

(2.5) T ld g = W k E Z er , T^G = {w \ w o e W k E « ler } V0G£, 

where W E '£ ler defined in ( II. ip is the phase space of the velocity fields of (E). From the Hodge 
decomposition, a complementary space of W E '£ ler in W k ' q (Q, M n ) is given by 

(2.6) (W k E : ier ) L = {Vh I h G H/ fc+1 ^,R)}. 
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Here the orthogonality is in the sense that 

/ w ■ Vhdx = 0, \/w G W k ' q ler , h e W k+1 ' q (Q, R). 
Jn 

It is clear that both W k ' q ler and (W^^) 1 are closed subspaces of W k ' q ({l, R n ) and 

w k ^w) = w k z er ®{w k / uler )\ 

In fact, given any X G W k,q (Q, M"), let 

(2.7) w = X - Vh 
where h is the solution of 

(2.8) Ah = V • X in V N h — X ■ N onffi, 

then obviously X = w + V/i with w G Wg' q ler and this verifies the direct sum. Locally 
near the identity map id, we will write Q as the graph of a smooth mapping from W E ' q ler to 
iy^Euier) 1 ' ano - thus ^ ^ s rigorously a smooth manifold. Let Bg(-) denote the ball of radius 5 
centered at in the corresponding Banach space. 

Proposition 2.1. There exists 5q > and a smooth mapping ^ : Bs (W E q ler ) — > W k ' q (fl, M. n ), 
such that (0) = id, £>tf (0) = I,V(w)-we (W k ' q ler ) L for any w G B 5o (W k ' q ler ) , and 

(id+Bs 4 (w k > q (n,R n ))^)ng c {*(«;) | ™ g B So (w k ' q ler )} c (id+B Mo (^(n,i n )))ne. 

Proof. Since <9f2 is a smooth compact hypersurface in M n , the distance function to d£l is 
smooth in a neighborhood of dQ. Let <i : M n — > R be a smooth function with compact 
support such that it coincides with this distance function in a neighborhood of dQ. Consider 
the mapping 

G : W k ' q {Q, R n ) -> y = {(/, #) G W^*- 1,9 ^, R) x Py fc ~H<9fi, R)\ f gdS = 0} 

Jon 

defined as 

G(0) = (det( J D0), (do^len--^- [ docjidS) 

where \dQ\ denotes the area of dil. Obviously, G is a smooth mapping and G\g = (1,0). 
Moreover, suppose </> G U and G(0) = (1,0) where U is a neighborhood of the identity map 
id in W k,q (Q, M ra ). Then is a diffeomorphism from Q to its image and 

\<f>(ft)\ — and d o = a = — ^— ■ / do<j)dS 

\<J^\ Jan 

imply that d o <p = a = 0. Otherwise a^0 would imply that (j>(fl) either strictly covers Q or 
is strictly contained in Q, either of which contradicts with the first identity above. Therefore 

ung = unG- 1 {(i,o)}. 

It is easy to compute 

VG{id)X = (V • X, X ■ N - I X-NdS)eY, VX G W k,q (Vl,~Ei n ). 
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From the standard theory of elliptic problems with Neumann boundary conditions, 



DG{id)Vh = (Ah, V N h - --j— / Ahdx) 



is an isomorphism from (W^ ler ) 1 ~ to Y . Therefore the proposition follows from the Implicit 
Function Theorem. □ 

Near any O G Q, the local coordinate map \I/ composed with the right translation, \l/(-)o0 o , 
gives a smooth local coordinate map near 0o- Therefore Q is a Banach submanifold with the 
model space W E '£ ler . It is well-known |EM70j that Q is not such a standard Lie group as one 
needs to be careful with the smoothness of the left translations. Let C(<f>\, 2 ) — 4>\ ° 02 and 
it is straightforward to verify 

Proposition 2.2. C G C l ([G n W^ fc+i ' 9 (f2, R™)) x and, /or any 2 G 0, iae rip/rf 

translation C(-,(p 2 ) e C°°(G). 

2.2. Euler equation as an ODE on TQ. It is well-known that the pressure p can be 
represented in terms of the velocity field v. In fact, by taking the inner product of v t and X 
on dVt and the divergence of the v t and using V • v = in Q and N ■ V = on dQ, we obtain 



(2.9) 



-Ap = Hl^d^dy) = Y^l^dividy = ti(Dv) 2 in Q 

\/ N p = -N ■ (v ■ V)v = -V v (v ■ N) + V V N -v = v U(v) on dQ 



where the symmetric operator II G C°°(fl,L(TdQ)) is the second fundamental form of dVl 
defines as IT(x)(t) = V t tX with r T G T x dQ being the tangential component of r. 

Based on the form of the pressure, we define the symmetric bounded bilinear mapping 
B : (W k ' q (tt, R™)) 2 -»■ W k ' q (tt, R n ) as 

£?(X 1) X 2 ) = V 7 

where 

' -A 7 = tr(DXiZ>X 2 ) = YIl )j=x diX{djX\ in Q 

Vtv7 = Xi ■ U(X 2 ) on dQ. 

The boundedness of B is clear from the standard elliptic theory. Note that here we do not 
assume V • Xi )2 = in Q or N ■ X± j2 on dQ in the definition of B. According to the Hodge 
decomposition given through (12. 8p and a similar calculation as in (12.91) . it holds that 

(2.10) DX 1 (X 2 ) + B{X U X 2 ) = (X 2 ■ V)X ± + B{X l: X 2 ) G W k / ulerl VX 1; X 2 G W k / uler . 

As a side remark, it indicates that, when embedded in L 2 (Q, R 71 ), B is the second fundamental 
form of Q at id which can be rigorously verified through a standard procedure. 
Define the mapping V : Q x (W k ' q (tt, R n )) 2 ->■ W k ' q (tt,R n ) as 

(2.11) V{<f>,X u X 2 ) =B{X 1 o<p~\ X 2 0-^0 0. 
We also define the projection Q : Q x W*>«(fi, R n ) -»■ T£ as 

(2. 12) Q(0, X) = (X o 0- 1 — V/i) o G 7^0 

where V/i for X o is defined in (12.71) . Obviously, V is symmetrically bilinear in X\ and 
X 2 . As in fl2TT0l . it holds 

(2.13) ((^(Xio^- 1 )) o0)(X 2 )+P(0,X 1; X 2 ) G T+Q = {w : Q ^ R" | wof 1 G W^}. 
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In fact, ( 12.1 Op also leads to that, when embedded in L 2 , V(<p, ■, ■) is the second fundamental 
form of Q at 0. Euler equation (E) and (I2.9P imply that the Euler equation (E) takes the 
form in the Lagrangian coordinates 

(2.14) u tt + V{u, ut, u t ) = 0, u(t)eG. 
Moreover, for any O £ G, we have 

(2.15) P(0o0 o ,X lO o ,X 2 o0 o ) =P(0,X 1 ,X 2 )o0 o , 
i.e. V is invariant under the right translation. Therefore, 

(2.16) u(t) o O is also a solution for any solution u(t) of (12.141) . 

The following proposition states that V is a smooth mapping and thus has no regularity loss 
which is far from trivial even though B is a bounded bilinear operator. To see this, one note 
that even the dependence of the term X\ o g W k,q on G W k,q is not smooth unless X\ 
belongs to a space of better regularity. The proof of the proposition is essentially a careful 
analysis of the commutator between the operations of B and the composition by G Q . 

Proposition 2.3. V : Q x (W k > q (tt,R n )) 2 W k ' q (Q,R n ) and Q : Q L(W k ' q (tt, R n )) are 
C°° and, for any m > 0, there exist Co, K > depending only on m and k such that 

(2.17) |(P ( ^) m 'P(0, Xi, X 2 )\ W k, q < C \D(j)\wk-i, q \X 1 \ W k, q \X2\ W k, q 

(2.18) \{V^) m Q{<P,X)\ W k, q < Co\D4>\K k -i„\X\ W H, q . 

Here {JD^) m V should be considered as a multilinear operator from T^Q to W k,q (Q). Since 
V is bilinear in X 12 , the bounds on the derivatives of V with respect to X 12 follow from 
(I2.17P for m = 0. Also, we note that \D<f)\ W k-i, q > 1 since k — 1 > ^ and detDcj) = 1. 

Proof. We will present only the proof for P as the one for Q follows through almost exactly 
the same (or even slightly simpler) procedure. Due to the invariance (12.151) of V under 
the right translation, we only need to show its smoothness near <p = id. As belongs to 
the manifold Q, the smoothness of V is equivalent to the smoothness of V(Q(w), X 1; X 2 J 
with respect to w G W E '^ ler and X i)2 G W k,q (Q, M. n ) for any smooth local coordinate map 
$ : B s {W^ ler ) -> Q. Moreover, to prove the Frechet smoothness of V, it suffices to show 
the Gateaux differentiability of V up to the m-th order for any m > 0, which would imply 
that Gateaux derivative T> m ~ 1 'P is continuous and thus it is also the (m — l)-th Frechet 
derivative. To show the Gateaux differentiability up to the m-th order, it suffices to prove 
the smoothness of 

P(0( Sl , . . . , s m ), X^st, . . . , s m ),X 2 ( Sl , s m )) G W k > q (Q, W 1 ) 

for any 0(s x , . . . , s m ) G Q C W k ' q {Q, R n ) and Xi )2 (si, ...,s m ) G W k ' q {Q, R n ) with smooth 
parameters (s±, . . . , s m ) G U C M m . We will show by induction 

(ML) Z = d Sl . ..d Sm V(4 > ,X 1 ,X 2 )\ Sl= ... =Sm=0 G W k ' q (Q,R n ) 

and obtain its bound in the form of (12.171) . 

The boundedness of the bilinear transformation B implies (12.171) and (ML) for m = 0. 
Assume (ML) and (12.171) hold for < m < mo and we will prove it for m = m^. Let 
P(0, Xi, X 2 ) = (V7) o 0, where 7, defined as in the definition of B, depends on sx, ■ ■ ■ , s mo 
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through 0, Xi, and X 2 . In the following, since it will be much easier to carry out some 
calculations in the Eulerian coordinates, let 

Xj = Xj o j = 1, 2; n = (d Si <p) o D Si =d Si +V n , i = 1, . . . ,m 

and 

Clearly, it is sufficient to show Z G W^' 9 , which will be achieved by studying its normal 
component on dQ, divergence, and curl. 

To start, for any vector field W G W k ~ 1+m °' q (£l, W 1 ), let Y(y) = or 
equivalently Dcf)(Y) = W o 0, then 

D0(9 ii y) = -D0 Si (y) + ((DW)o0)0 a4 d s Y eW k -^(Q,R n ). 

Differentiating the above identity one more time implies 

D(j)(d SjS Y) = -Dcj) s .{d s Y)-D<l> Si {d s Y)-D^^ 

and thus the smoothness of and W yield 

9 Sj /^ H? (0,l n ) ifm >2. 

Repeating this procedure we obtain d Sii ,„ 8i Y G W /fc_1,l3 (f2, IR n ) inductively. Changing from 
the Eulerian coordinates to the Lagrangian coordinates, we have 



(V W Z -D S1 ... D Smo (D 2 j(W))) o<j ) = V Y Z-d Sl ... d Smo Vy ((V 7 ) o 

Applying the induction assumption to the last term above and using the commutator formula 
[d Si , Vy] = Vd Si y to move the Vy to the outside to produce VyZ, we obtain 

(2.19) (V W Z -D S1 ... D^D^W))) o<p = V Y Z-d Sl ... d Sm( V Y ((V 7 ) ° 0) G W k ~^. 

Taking W = e±, . . . , e n of the standard basis of lR n , (12.191) and the definition of V imply 
the curl V x Z contains only the commutators terms and thus satisfies 

V x Z G W* _1 » a . 

Similarly the divergence satisfies 

V-Z + D S1 ... D Smo (tv(DXi) (DX 2 )) G W k ~^. 

Expanding the above last term using the product rule on D s ., it consists of terms in the 
form of 

(D Sn ...D Sim d h X[ 2 ){D Sn ...D Sjmo _ m di 2 X l 2 1 ), {i'i,...,i m , j h . . .j mo - m } = {l,...,m }. 

Move and di 2 to the outside in the same fashion as in the derivation of (I2.19p (replacing 
W by e\ x 2 and V7 by X 1 ^ ) and using the smoothness of X\ % in si, . . . , s mo , it is easy to 
obtain D Sl . . . D Smo (tr(Mi)(M 2 )) G W*- 1 * and thus 

V • Z G W^ -1 ' 9 . 

Finally, using 

D Si N = Ii( Ti ), D Sj D Si N = (V Tj II)(7i) + IT(0 s , Si o 0- 1 ), . . . 
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the induction assumption, and the assumption iV • V7 = on dfl in the definition of V it is 
straight forward to obtain 

Z ■ N = N ■ D S1 . . . D Smo V 1 G W k -«> q (dn, R) 

in the same fashion. Therefore, the standard estimates in elliptic theory implies that (ML) 
holds for m = m . Moreover, inequality (12.1 T[) follows from the observation that the compo- 
sition by <p or cf)" 1 only produces terms like \<P\ l w k m ^ ne estimates of the W k,q norms and 
thus the proof of the proposition is complete. □ 

Proposition 12.31 provides the key element for us to prove that (12.141) is a smooth second 
order ODE on the infinite dimensional configuration manifold Q. 

Proposition 2.4. For any uq G Q and Wq G T U0 Q , the initial value problem of the Eu- 
ler equation (12.141) has a unique solution (u(t),ut(t)) G TQ, locally in time, depending on 
(t, Uo, Wo) smoothly. 

Proof. Due to the right translation of (12.141) given in (12.161) . we may assume u$ belongs to a 
small neighborhood of id. From Proposition 12. 1[ 

(w,Vh) -> *(«;) + V/i 

is a local diffeomorphism from W E '^ er x (W^^) 1 - to W k ' q (Q, M. n ). Taking the ^(w) compo- 
nent, we obtain a smooth $ : id + B s (W k ' q {Q, R n )) ->■ Q such that 

(2.20) w-$(u) G (W^ ler ) x , Vueid + B 5 (W k ' q (n,m. n )) 
and 

$(«) = u, V« G gn (id+ B s (W k ' q {Q,R n ))y 
Consider a modification of (12. 14|) 

(2.21) utt + V(§{u) 1 u u V<b{u)u t ) =0, u G B 5 (W k ' q (n,R n )) u t G B s (W k ' q (Q, W 1 )) . 

From the smoothness of $ and V, equation (1 2 . 2 1 [) is a smooth ODE defined on an open 
subset of an infinite dimensional Banach space (W k,q (Q, M n )) 2 and thus is locally well-posed 
with smooth dependence on the initial value. Moreover, any solution u(t) of (12.211) satisfying 
u(t) G Q for all t also solves (I2.14p . Therefore, to complete the proof, we only need to show 
that, if the initial data is given on TQ, then the solution of (12.21 j) also stays on TQ, i. e. 
u(t) G Q and u t (t) G T u ^ t )Q. Let v = u t o ($(tt)) . Equation (I2.21|) yields 

iH o $(u) + ((Dv) o (V${u)ut) = u tt = -V($(u),u t , V$(u)ut) 

and thus (12. 13[) implies v t o g T$( u )^ or equivalently v t G = W E ' q ler . Therefore 
v(t) G W E ,q ler , and thus u t = v o $(m) G X$( u )(/, for all t follows from the initial assumption. 
Consequently 

(u - $(u)) t = (J - X?$(m)H = 

where we also used that (I2.20j) implies VQ(u)X = X for any X G T<&/ U )G- Therefore 
u = $(tt) G Q and the proof completes. □ 

According to (12.161) . the second order ODE (12.141) defined on the Lie group Q is invariant 
under the right translation. The standard procedure of taking v = Ut o u^ 1 reduces it to a 
first order equation on the corresponding Lie algebra W E ' q ler = T id Q, which turns out to be 
the usual form (E) of the Euler equation in the Eulerian coordinates. However, according 
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to Proposition 12. 2[ the composition as the multiplication operation on the group Q is not 
smooth, this procedure induces the loss of one order of spatial derivative and make (E) into 
a PDE, i. e. the right side of (E) does not define a smooth vector field on W E ^ ler . 

2.3. Euler equation near Vq as a non-autonomous ODE on TQ. Even though Euler 
equation is equivalent to an infinite dimensional ODE on TQ, the Lagrangian frame work 
also brings two complications: 

• Q is not flat, which means that we may have to carry out the analysis in local 
coordinates on Q and 

• the steady velocity field vo(x) of the Euler equation (E) corresponds to a dynamic 
solution (u (t,y),u t(t,y) = v (u(t,y))) of (|2.14|). 

It is natural to look for ways to take the advantage of the group structure of Q to reduce 
(I2.14p . localized near uo(t), to a non-autonomous second order ODE defined in local coor- 
dinate neighborhood of id. Based on the comments at the end of Subsection 12.21 taking 
w = u o Uq 1 for u(t) near uo(t) would result in loss of regularity, which can also be seen ex- 
plicitly through the simple calculation u t = w t o uq + (Dw o uq)(vo o uq) leading to w t £ W k ' q 
due to the presence of Dw. Instead, for any solution (u(t),u t (t) = v(t) o u(t)) G TQ of the 
Euler equation (I2.14p with u(t) close to u (t), let 

(2.22) « = *(*, to) 4 «,,(*) <,*(«;(*)), w(t) G S fc (Wj$J, 

where $ is given in Proposition 12.11 which also implies <&(£,-), for any t G R, is a local 
diffeomorphism from W E ^ ler to Q. Therefore, the linearizations 

X = V$(t, w)X = (Du ° *H)wwz, X e W^ ler 

X = (V^w^X = (V^iw^dDuo)- 1 o *(w))X, X G T U Q 

are isomorphisms from between W E £ ler and T U Q. Substitute (12.221) in to (12.141) . one may 
compute 

(2.23) u t = $ t (t, w) + V$(t, w)w t = u ot o + (Duq o ^(w))V^(w)w t 

(2.24) v = u t o u- 1 = v Q + ((Duo) o Uq 1 ) ((V^(w)wt) o ^(w)' 1 o Uq 1 ) 

and 

utt =$tt(t, w) + 2V$ t (t, w)wt + V$(t, w)w tt + P 2 $(t, w)(w t , w t ) 

=u oa o $?(w) + 2 (Du ot o ty(w)) W(w)w t + (D 2 u o ty(w)) (W(w)w u W(w)w t ) 
+ (Dm o$H) (V 2 y(w)(w t ,w t )+VV (w)w tt ) . 

Therefore, for u(t) close to u (t), the Euler equation (I2.14p is rewritten as 

(2.25) w a + T(t, w, w t ) =0, we B 5o (W k E « er ), v> t G W^, 

where, for w G B So (W^ ler ) and X G W^ ler , the term F(t,w,X) is arranged into the linear 
and quadratic parts in X 

(2.26) F(t,w,X)=A(t,w)X + B(t,w)(X,X) G W** er 
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with the linear and bilinear (in X) operators A(t,w) and B(t,w) are given by 
A(t, w)X =2 (Z>$(i, w)) ^ (V$ t (t, w)X + V{u, v o u, X)) 

(2.27) =2 (£>tf ((Duq)- 1 ° # W) ((Dv o u)X + V{u, v o w, X)) 

5(t,w)(X,X) =(P$(t,w)) _1 (P 2 $(t,w;)(X,X) +P(m,X,X)) 

= (P^(w)) _1 (( j Dm )- 1 o V(w)) ((Du o ^(w))X? 2 ^(w)(X,X) 

(2.28) + (D 2 u o *(«;)) (D*(w)X, Wf{w)X) + V(u, X, X) 
where 

U = $(t, iy) = u o *(iy) X = V$(t, w)X = (Du o ty(w))W(w)X. 

In the above calculation, the invariance of V under the right translation (I2.15P and equa- 
tion (I2.14p were used in handling both u a and Uq U - Here even though the linear operator 
w)) acts only on the subspace T U Q, the terms A(t,w)X and B(t,w)(X, X) are 
well-defined and thus T(t,w,X) G W^ ler . In fact, (I2.13P implies 

£>$t(t, w)X + V(u, v o u , X) = (Dv o u)X + V(u, v o U) X) G T U Q 

and thus A(t, w) is well-defined. To see that B(t,w) is well-defined, we note the second 
linearization of $ along a line w + sX, s G R, at s = 0, is given by 

m ss = V 2 <S>(t,w)(X,X) = (Du o^(w))V 2 ^(w)(X 1 X) + (D 2 u o^(w))(V^(w)X,V^(w)X). 
Let 

F^.or^io^e W k E * er , 

then 

M ss = y s o M + J DF(X). 
Since K s G Wj$ er and (EUJ) implies 

DF(X)+P( M ,X,X) e^, 

we obtain 

P 2 $(t, w)(X, X) + +P(w, X, X) = u ss + V(u, u s , u s ) G T U Q 

and thus 

B(t,w)(X,X) = (V$(t,w))- 1 (V 2 $(t,w)(X,X)+V(u,X,X)) E W^ ler 
is well-defined. 

Remark. An alternative way to rewrite the Euler 's equation to derive the above form is to 
follow the Lagrangian variational principle. One may first express the action J J n ^-dydt, 
defined on TQ, using ( 12.22ft and ( 12.231) . The Euler's equation in terms of w follows from the 
the variation of the action. 

Recall we assumed in (Al) that Vq G W k+r,q with r > 4 and fi > /xq > is a constant 
fixed before (1231). 
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Lemma 2.5. The nonlinear mapping J 7 satisfies 

7 G (R x B So (W k E ZJ x W k E « ler , W k Z er ) , Ht, w, 0) = 0. 

Moreover, there exist K > depending only on r and k and C > depending only on 
r,n,k,q,v , such that, fort 6R, «j,I6 ^Euier an< ^ \ w \w k >i < ^0; 1,76 ^ a,ue 

i /?<v .11 / \ < rv-^'*! 

1 ^ ' jl C-2(B, (^ e J,L(^ ier ®^ er ,H/^)) " 

i/ere i/ie highest order derivative is in the Gateaux sense which is sufficient to yield the C r ~ 3,1 
bounds. 

Proof. The smoothness of T follows directly from its expression f !2.26|) - ( I2.28|) and Propo- 
sitions [2J] - [231 The property ^(t, w,0) = follows directly from (I2.26j) which is actually 
a consequence of the right translation invariance of the Euler equation. To demonstrate the 
latter, we notice that (u o <fi, u ot o 0) is a solution of (I2.14p for any G Q. Then (12.221) 
implies that, for any w G Bs (W E ^ ler ), (w,0) is a time independent solution of (I2.25P and 
thus w, 0) = 0. The derivation of the estimates is tedious, but straightforwardly from 
Proposition O-ESl (12~T7D . and fl2T26l) - (I2T281) . 

□ 

Remark. ^4s proved in Proposition \2.3[ (I2.14p a smooth infinite dimensional ODE. How- 
ever, the local coordinate systems based on the composition would always cause loss of deriva- 
tives due to Proposition \2.2\ Here our local coordinate mapping <3>(t, •) allows us to obtain 
some limited smoothness due to assumption (Al) of the extra regularity of vq. 

2.4. Linear exponential dichotomy in Lagrangian coordinates. Since Fit, w,0) = 
for any small w G T id Q, we have V w F{t, w, 0) = 0. We can rewrite (I2.25P as 

(2.29) z t = A (t)z + F(t, z), z = (z u z 2 ) T G B So (W^ ler ) x W k / uler 



where 



Mt) ~ (o -A{t, 0)) ' Z) ~ \ A(t, 0)z 2 z u z 2 ) 



From (I2.27p . the explicit form of A(t, 0) is given by 

A{t, 0)X = 2{Du {t)Y l [(Dv o u {t))Du {t)X + V(u , v o u {t), Du {t)X) 
and Lemma [2.51 implies that there exists C > independent of So such that 

(2.30) F(t,0) = = P M VlFit, < CS* 

where again the highest order derivative is in the Gateaux sense which is sufficient to yield 
the C r_3,1 bounds. The linearization of (12.251) takes the form of 

(2.31) w tt + V x 7(t, 0, 0)w t = 0^z t = A (t)z 

whose well-posedness is guaranteed by Lemma [2.51 Let T(t,t ) be the solution operator of 
(I2.3ip with initial time to and terminal time t. 
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On the one hand, for w G W^ ler , z = (w, 0) T is a solution of ( I2.3ip . On the other hand, 
linearizing ( 12.241) at the steady solution z = (wo,0) T one may compute z = (w(t),w t (t)) T is 
a solution of the linearization of ( 12. 29ft at z = (wq, 0) t , where 

w t (t) = (VV{w )) 1 ({Duoit))- 1 o *{w )) (v(t) o uo(t) o *(«,„)) 

and f = is a solution of ( 11.2ft . In particular, taking wq = 0, we have that z = 
(w(t),w t (t)) T is a solution solution of ( 12.311) where 

(2.32) w t (t) = (Du (t)) 1 (v(t) o Uo (t)) . 

This correspondence between the linearized solutions of (II. 2p and (I2.3ip and assumption 
(A2) would yield the (t-dependent) exponential dichotomy (W^ ler ) 2 = Y u (t)®Y cs (t) of (EOT]) 
such that T(t, to)Y UtCS (to) = Y UfiS {t) and the exponential decay rate of T(t, t ) as t — > — oo (or 
+oo) in Y u {t) (or Y cs (t)) is bounded roughly by \ u (or A cs ). To define Y u (t), it is natural from 
(I2.32p that, for (w,w t ) G Y u (t), w t takes the form of (Du (t)) 1 (e tL v o u (t)) with v G X u , 
where L (as well as L u>cs ) is the linear operator defined in the linearized Euler equation (II. 2ft . 
Also the decay of w as t — > — oo requires it take the form of w(t) = f_ Wtdt'. Therefore, 
for any t el, let 

K(*) = {( / {Du Q {T)Y\{e TL v)ou {T))dT, {Du Q {t)Y\e tL vou {t))) T \v&X u }. 

J — oo 

The convergence of the above infinite integral follows directly from assumptions (A2) and 
(A3) with K sufficient large depending only on k. Similarly, for (w,w t ) G Y cs (t), w t takes 
the form of (Duo{t)) (e tL v o uo(t)) with v G X cs , and w should take the form of 

w + [ (Du (T)y 1 ({e tL v)ou {T))dT. 
Jo 

However, w G W E '£ ler is arbitrary and thus we can absorb the integral term into w . Define 
Y cs (t) = {(w, {Du {t)Y l {e tL vou {t))) 7 \weW^ ler , v G X cs }. 

Lemma 2.6. It holds (W^J 2 = Y cs (t) © Y u (t) . Moreover, let P u , cs (t) G L((W^ ler ) 2 ) ^ 
the projections associate to this decomposition and 

T UjCS (t,t ) = T(t,t )\ YuiCa (t ). 

Then for any t, t < 0, we have T UtCS (t, to)Y UtCS (t ) = Y UiCS (t) and there exist constants K > 
depending only on k and Co > 1 depending only on k, n, q, Vo such that if \ u > Kfi, we have 

\T cs (t,t )\ < C e A -('-'°)-^ , V > t > t 
\T u {t,t )\ < c e (A "-^ )( '-' o) -^'°, V t < t < 0. 

Remark. This lemma shows that unlike the traditional exponential dichotomy, the norms of 
the projections in the invariant splitting here is not uniformly bounded in t and may approach 
oo as t — > — oo. This means that the angles between the unstable and center-stable subspaces 
Y UjCS {t) of (I2.3ip may not have a uniform positive lower bound as t — )■ — oo. 
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Proof. We first show the invariance of Y cs>u (t) under T(t,t ). In fact, for any 

Y u (t )3 z=(w, Wl ) T = ( f° (Du (T))~\(e TL v)ou (T))dT, (Duo^))* 1 {e toL v o u (t ))) T , 

J — oo 

where v G X u , let 

v(t) = e tL v G X u , w= (Du {r)y 1 ({e TL v)ou {T))dr. 

J —oo 

Clearly, 

Y u (t) 3 z(t) =(w(t) lWl (t)) T 4 f f {Du {r)Y\{e TL v)o Uo {T))dT, (Du^t))' 1 (e tL v o u (t)) 

J —oo 

= (^,0) T +(^ ( J Dn (r))" 1 (t;(r)o Mo (r))dr, (^(tjj^^o^fi)))' 

is the solution of (12.311) with initial data (w,v) at t — 0. Moreover, it satisfies z(to) = -2 and 
thus Yu(t) 9 = T(t,t )z which implies the invariance of Y u {t). 

The above arguments also leads to the decay estimate of T u (t, t ). In fact, one may compute 



Wl (t) = {Du (t)) "^-^((DuMw!) o K(to))- 1 )) ou Q (t) 
Since v G X u , we obtain from (12. 3p and assumptions (A2) and (A3) 

I^iWIwm < Ce^^^-^Vil^,^ Vt < t < 0. 

Finally from w(t) = f_ Wi(r)dT, we obtain the estimate for T u (t,to). The proof of the 
invariance of Y cs (t) and the estimate for T cs (t,t ) are similar. 

To prove the direct sum and obtain the bounds on the projection operators, let P c ° s u G 

^(^Euier) b e the projections given by the decomposition W^ ler = X cs © X u assumed in 
hypothesis (A2). Given any z = (w,wi) T G (W E ^ ler ) 2 and t < 0, let 

v C s,u(t) = ^cV^ L (( Du o(t)^i) °u (t)- 1 ) G X CS)U 
w (t) =w- f (Du {r)Y l ({e TL v u {t)) ou (r))dr 

J — oo 

From (12. 3p and assumptions (A2) and (A3), we have 

|e T V(*)lw*.« < Ce^^-^Vil^,,, V r < t. 

Let 

T 



' {Du (t)) \{e TL v u {t))ou Q {r))dr 1 (Du (t)) 1 (e tL v u (t) o u (t))) eY u (t) 

«/ — oo 

Zcs =(w (t), {Du {t)Y\e tL v cs {t) ou {t)Y) T G Y cs (t). 

Obvious z = z u + z cs and this splitting is unique. Therefore {W^ ler f = Y cs (t) © Y u (t) and 
Pcs,u{t)z = z CSjU . It is straight forward to first obtain the estimates on P u (t) based on the 
above inequalities and the bound on P cs {t) = I — P u (t) also follows. □ 
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Let 

Fu,cs(ty Pu,cs(t)F'(t,Z). 

Then Lemma 12.61 and (12.301) imply that there exist K > and G\ > such that for any 
t < 0, it holds 

(2.33) F cs u (t,0) = = T) Z F CS u (t, 0), \T) F ucs (t, z)\ rr _ i(R (w h, q x 2 w k ' q i — ^*i e ^ > 

where again the highest order derivative is in the Gateaux sense which is sufficient to yield 
the C"* -3 ' 1 bounds. 

2.5. Integral unstable manifolds of (I2.25p . We will follow the Lyapunov-Perron integral 
equation method to construct the integral unstable manifolds. In the standard construction 
of invariant manifolds, where the bounds on the invariant splitting and the nonlinear terms 
are uniform in t, small Lipschitz constant of F csu is sufficient in the construction of local 
invariant manifolds. As we do not have t-uniform estimates here, we repeatedly used the 
property that F(t, z) = 0(\z\ 2 ) to yield an extra \z(t)\ = 0(e xt ) with At < 0. This quadratic 
nature of F combined with the exponential gap condition (I2.36P allows us to complete the 
proof of Proposition 12.71 in the below. However, similar construction would not work for the 
construction of the center-stable manifold since solutions z(t) on the center-stable manifold 
do not satisfy z(t) — > as t — > ±oo. 

For A G (A cs , X u ) and <5i < Sq to be determined, let 



r = {z= (z u ,z cs ) G C ((-oc,0],B Sl (W k E ^ ler r) | \z\ x < 5,} 
where z U)CS (t) G Y UiCS (t) and 

\z\ x = sup e~ xt \z(t)\ W k, q . 

This r is the set of functions with the desired backward in time decay expected to be satisfied 
by the solutions on the unstable manifolds. For any z G r and 

5x 



(2.34) z u0 G y„(0), 



2C n 



where C is given in Lemma [276| define £(-,z u0 ), where z = (z u , z cs ) = C(z,z u0 ) is defined 
as, for any t < 0, 



(2.35) 




T u (t, 0)z u0 + f*T u (t,T)F u (r, z(r))dT 
~- flooTcsit, t)F cs (t, z(r))dr. 



Proposition 2.7. There exists K > depending only on r and k such that if A and Si 
satisfy 

(2.36) A G (X cs + Kfi, X u - Kfx) 

(2.37) CoC&i-r-Xr- + t 1 r) < \ 

A - A cs X u - Kfx- y 2 

where Cq and C\ are from Lemma lKUi and (12.331) . respectively, thenC(-,z u o) is a contraction 
on r with the Lipschitz constant | for any z uQ satisfying ( I2.34[) . Moreover \£\ C r-3,i < C for 
some C > depending only on r, k, q, n, and Vq. 

Remark. Assumption (A3) with a reasonably large K , depending only on r and k, guarantees 
the existence of X and 5\ satisfying the above inequalities. 
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Remark. It is standard in the invariant manifold theory (see, for example, |CL88j ) to prove 
that z G r solves (12.291) with z u (0) = z u0 if and only if, z is the fixed point of C(-,z u0 ). 
Therefore, Proposition \2. 7| shows that, for any given z u0 satisfying (12. 34ft . there exists a 
unique solution of (I2.29p satisfying the exponential decay as t — > — oo with the decay rate at 
least A. 

Proof. In the proof the generic constant K may change from line to line, but always depends 
only on k and r. From the definition of C, Lemma [2761 assumptions (A3), (12.361) . and (I2.34p . 
and the second order Taylor expansion of F based on (I2.33p . (instead of the usual small 
Lipschitz estimates of F), we obtain for t < 



e- xt \~zu(t)\ w ^ <C \z u0 \ wk , g + y° ic C ie - A '+( A «-^)('-)-^+ 2A ^r| 

<C \z u0 \ wk , g + \c Cl\z\\J e (A«-^-A)(t-r) dT 



*-n I I CqC x 5i 

<c KokM + 2(Au _^_ A) |^ 



and 



-Ati~ 



CqCi8i\z\\ 
2(A - \ cs ) ' 



,(t)\w>*< [ l-CoC^^-^^drlzllK 

J -oo " 

Therefore (I2.37P implies that 

1 3 

(2.38) \C(z, z u0 )\x < C \z u0 \ W k, q + -\z\ x < ~5i 

and thus £(•, z u o) maps r into itself. 

To prove £(•, z u q) is a contraction on J 1 , we note that for any z 1,2 G -T, (I2.33P implies 

|F U)CS (M 2 (0) - F U!CS {t,z\t))\ < CAe^^z 2 - z l \ x , Vt < 
Therefore, we have 



l(t)\ <GA*i ^ e^-^-^dr\z 2 - z\ < °f ^ - ^ 
Jt K — Kfj, — X 

e- Xt \zUt) - <C Q C& f e^-^dr\z 2 - z\ < 



and thus (12.371) implies that C(-,z u0 ) is a contraction. 

Since C is linear in z u q, we only need to prove its smoothness in z G r. For any z G -T, 
formally the linearization of £ is given by 

T> Z C{Z, Z u q)Z\ —Z\— [Zlu, Z\ cs ) 

where for t < 0, 

ziu(t) = jQT u (t,r)V z F u (r, z(T))zi(r)dT 
zics{t) = J_! oo T cs (t,r)r> 2 F cs (r,2;(r))2; 1 (r)c/r. 

The same procedure as in the above shows that |I> 2 £|a < \- To show it is indeed the 
derivative of £, take G -T, let 

Z\,2 = {Zl,2u, Z\,2cs) = £(Zi t 2, Z m q), Z = (z u , Z cs ) = V z C(zi, Z u q){z2 — Z\) . 



(2.39) 
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It is straight forward to compute from (I2.33j) and Lemma [2.61 that for t < 0, 
e- xt \z 2u {t) -z lu {t) -z u (t)\ 



-e- xt 



J T u (t,T)^F u (T,z 2 (T))-F u (T,z l (T))-VF u (T,z 1 (T))(z 2 (r)-z l (r))yT\ 



<CqC\ \Z2 — Z\ 



1 2 / , -A/-I 



n n lr 2 _ r 1 ! 2 

t+{X u -K l M)(t-T)-K f J,T + 2XT ^ < MDWlf z lA 

~~ A u — fT/i — A 

The estimate for the center-stable component is very similar and this proves that T> Z C is 
indeed the derivative of C. 

Finally, we will show that D Z L is Lipschitz. In fact, let 

Z\ % 2 = {Z\2ui Z\2cs) — Dz£(Zi t 2, Z u0 )z. 

Then for any t < 0, 

e- xt \z 2cs (t) -z lcs {t)\ = e~ xt \ I T cs (t,T)(V z F cs (T,Z2(r))-V z F(T, Zl (T)))z(T)dT\ 

J — oo 

<QAh - Zl \ x \z\x f e^-^-^- K ^dr < y^\z 2 ~ zi\x\z\x. 



oo 



The estimates for the unstable component is the same and the proof of the higher order 
smoothness is similar. □ 

From the Contraction Mapping Theorem and Proposition 12. 71 the mapping £ has a unique 
fixed point 

which is C r ~ 3,1 in z u q and satisfies -2* u (0,z u o) = z u q. Moreover, (12.381) implies 
(2.40) \z*(u )\x < 2C \z u0 \ W k, q . 

Like in the standard invariant manifold theory, these are solutions on the invariant integral 
unstable manifold of the non-autonomous system (12.251) . Define 

^-l(^mo) = £*es(0, Z u q) G X cs 

for all z u o satisfying (12.341) . then the C r_3,1 graph 

Wl 4 graph(h L ) 

defines the slice of the unstable integral manifold of (I2.25P for i — 0. Obviously the unique- 
ness of the fixed point implies that z*(-, 0) = and thus h L (0) = and G W^. Differenti- 
ating the fixed point equation we obtain 

T^z u0 z *{z u q) = D ZuQ C(z*(z u o), z u o) + T> z C(z*(z u q), z u o)'D Zu0 z*(z u o) 

Clearly, (I2.39P implies that V z C(0, z u0 ) = and thus 

V ZuO z*(-,0) = V Zu0 £(0,0) = (T u (;0),0) 

which does not have the center-stable component. Therefore, we obtain that 

MO) = o 

which means that, at 0, the tangent space of the unstable integral manifold Wl 

T Wl = K(0) = {(U(v), v) T | v G X u }, 
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where 

(2.41) U(v) 4 I (L>« (r))~ 1 ((e^) ou (r))dr, \U\ L{Xu>wk , ) < oo. 

J— oo 

Here the boundedness of U follows from ( 12. 31) . (I2.36j) and assumptions (A2) and (A3). 

2.6. Unstable manifold in the Eulerian coordinates. From the unstable integral man- 
ifold (at t = 0) W1 constructed in the Lagrangian coordinates and the corresponding rela- 
tionship given in (12.221) and (I2.24p . we obtain the C r ~ 3,1 invariant unstable manifold in the 
Eulerian coordinates 

W u = {v = v + (V^f(w)w 1 ) o ^(w)- 1 | (w, wx) G Wl}. 

The above expression was derived by substituting t = into (12.24)) . 

Lemma 2.8. There exist K > (depending only on r and k), 82, C > depending only on 
r, k, q, n, and vq such that 

(1) There exists H G C r - 3A (B 3 s 2 (X u ), X cs ) satisfying \H\ C r- 3 ,i < C, H(0) = 0, DH(0) = 
0, and 

{v +v 1 +H(v 1 ) I vx G B S2 (X U )} C W u n(v +B 2S2 (W^ ler )) C {v +Vx+H{vx) I v x G B 3 s 2 (X u )}. 

(2) For any v# G (v + B 2 s 2 (W #^ er .)) ^ W u the solution v{t) of the Euler equation (E) 
with the initial value v (0) = v# satisfies 

v(t) G W u , \v(t) - v \ W k, q < C\v # - v \ wk , q e ix - K ^, V t < 0. 

Proof. (1) For any G X u with 

. I < S 1 

2c {i + \u\ L(x ^ Kzj y 

let z u0 = (Uvx,Vx) G Y u (0), where U is defined in ( 12.4 ip and it implies z u0 satisfies ( 12.34)) . 
and 

G{vx) = (V^{w)w x ) o ^(w)- 1 where (w, w x ) = 2,(0, ^ u0 ) = z u o + h L {z u0 ) G 

Clearly the definition of W% and W u and the properties of hi, imply W u = {vq + G(v x )} and 

G(0) = 0, DG{0)vx = w 1 = vx, 

and G G C"* -3 ' 1 with bounds depending only on r,k,q,n, and Vq. Therefore, the existence 
and properties of H follows from the Implicit Function Theorem immediately. 
(2) For any 



V # = ^0 + O G («o + S Ma (^er)) n W ^ With ( W #> ">!#) G W , 



L ■ 



let f(t) be the solution of (E) with v(0) = v# and z(t) = (w(t),wt(t)) the solution of (12.251) 
with the initial value z(0) = (w#,wx#)- Since z(0) = (w#,Wx#) G W% and thus z G r and 
023QD implies 

|w(t)|w*>9 + < e *kU < 2C e * \z u0 \w h ,i, t < 0. 

Therefore (I2.24p implies the desired decay estimate. 

\v(t) - v \w^ < Ce (x ~ KlM)t ->■ as t -> -00. 



UNSTABLE MANIFOLDS 21 

To see the local invariance of W u , fixed T > 0, for t G (— oo, T], let 

w(t) = tf" 1 (uo(t Q ) oV(w(t + t )) o uo(to)" 1 ) , t < 0. 

Since the right translation invariance and (12.221) imply 

u(t) = Uo(t) o (w(t)) = u (t + t ) o #(w(t + t )) o uo(to)^ 1 = u(t + t ) o M (t ) _1 

is a solution of (I2.14p . we have (w,w t ) is a solution of (12.251) which clearly corresponds to 
the solution v(t) = v(t + t ) of the Euler equation (E). Moreover, Proposition 12.11 implies 

\w(t)\ W k,q + \w t (t)\ W k, g < C(\w(t + t Q )\ W h,q + \Wt(t + t )\ W k,q) < Ce ) | Z u0 \ W k,q 

< Ce XT {l + \U\ L{x ^ Kz J \v# - v \ wKq e xt 
<2^e^(l + |^| w ,, u ) e ^. 

By choosing 

5,< h 



we obtain that the solution z(t) = (w(t),w t (t)) £ r. Therefore v(to) = v(0) £ W u which 
implies the invariance. □ 

The property DH(0) = immediately implies, as expected, the tangent space at the 
steady state vq is given by 

T V0 W U = X u 
and the proof of Theorem 11.11 is complete. 

3. Two-dimensional Euler equations 

In this and the next sections, we will illustrate how assumptions (Al) - (A3) can be 
satisfied for certain steady states of the Euler equation (E). In this section, we consider the 
case of fl = S 1 x (—yo,Vo), that is, 2ti— periodic in x and with rigid walls on {y = ±y }- 

Let v = (v i, t>2) T : — > ffi 2 satisfy V • v = in Q and v ■ N = on dQ. On the one hand, 

let 

1 f 

u = d x v 2 — dyVi, s = I v\dxdy 



|0| „ 

be the curl and the average horizontal momentum, respectively. We note that s is an invariant 
of (E) due to the translation symmetry in x. On the other hand, v is uniquely determined 
by lo and s through 

(3.1) v = JVA- 1 cu + se 1 , J = ^ , e 1 = (l,0) T 

where A -1 is the inverse of the Laplacian with zero Dirichlet boundary condition. It is clear 
that 

— \v\ W k, q < \w\ W k-i, g + \s\ < C\v\ W k, q , k>l, q > 1 
for some C > 0. In the (u, s) representation, (E) takes the form 

(3.2) u t + v -Vu = 0, s t = 
where v is considered as determined by (u,s) by (13. ip . 
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Remark. Due to the nontrivial first cohomology group of Q, the vorticity alone does not 
determine a vector field in W E '^ ler and thus the average horizontal momentum has to be 
included in the reformulation of the problem. If one considers Q = T 2 , the 2D torus, then 
both momentum invariants each of which corresponds to a nontrivial element in the first 
cohomology groups should be included. 

Suppose v G W k+4,q , k > 1 + |, is a steady state of (E), which corresponds to (o>o, so) 
has Lyapunov exponent > (both forward and backward in time). We linearize (13.21) at 
(ujq, so) to obtain 

(3 ' 3) W«-"l ° )~\ )- L "{s) +L ^s 

Assume that there exists an unstable eigenvalue Ao with ReAo > {k — 1) /i of the lin- 
earized Euler operator L (defined in (II. 2p ) on L q , and let v G L q be the eigenfunction with 
corresponding u and s. Then obviously s = and 

X u + v ■ \/tu = —v ■ Vwo- 

An integration of above along the steady trajectory X (s) yields 



/"OO 

u= / e-H-Vu (X (s))(ls. 
Jo 



By the standard bootstrap argument and the assumption ReAo > (k — 1) fio, we get uj G 
W k ~^ q and v G W 7 ^ 9 . 

For any A_ G ((A; — l)/io, Re A ) which does not equal the real part of any unstable eigen- 
value , let 

X cs = {{u,s) T | u G W k ~ l > q | Hmsup-logle^+^wl^-M < A_} 
which is clearly a invariant subspace of e Lo+Ll . Let 

a cs = a(e Lo+Ll \x cs ), o u = a(e Lo+Ll \ wk -i, qxR )\a CS} A+ = log(inf{|A| | A G a u }) > A_. 

As the groups of bounded operators e *( io+Ll ) and e tL are conjugate through (13. ip . we also 
have the invariance of X CS)U under e tL . Since e tL °u = uo« () (t)" 1 ) inequality (12. 3 p implies 
\ etL °\L(w k - 1 -i) ^ Ce^ -1 ^'*' for any /i > /i and some C > depending on /i. In particular, 
f is divergence free, yields that e tL ° is a group of isometries on any L q space. Since Lx is a 
compact operator acting on (co, s) T , e*( L ° +il ) is a compact perturbation to e tL ° in the space 

• A + > A_ and a u is an isolated compact subset of o"(e Lo+il ). Let X u be the eigenspace 
of e Lo+Ll corresponding to a u , and 

X cs>u = {veW^ ler \(uj,s)eX CS!U }. 

• X u and X u are finite dimensional and 

• (A2) is satisfied for any X cs and A„ with A_ < A cs < \ u < A + . 

Assumption (A3) depends on the Lyapunov exponents of vq. In particular, if vq is a 
linearly unstable shear flow, fio = and (A3) is also satisfied. An example is vq = (sin/3?/, 0). 

By |L03t Theorem 1.2], vq is linearly unstable when (5 > 1 and (^j^j < f3 2 — 1. 
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Remark. Consider rotating flows v = U (r) eg in an annulus Q = {a < r < b}. Then by 
similar arguments as above, assumptions (A2)-(A3) are satisfied as long as v is linearly 
unstable. 

4. Three-dimensional Euler equations 

In this section, we construct examples of 3D unstable steady flows for which Theorem 11.11 
can be applied to get unstable (stable) manifolds. Consider Q = T 3 to be a 3D torus with 
periods L x ,L y and L z in x,y and z variables. For any profile U (y, z) , the 3D shear flow 
uq = (U (y, z) , 0, 0) is a steady solution of 3D Euler equation. We construct unstable 3D 
shears satisfying assumptions (Al)-(A3) in several steps. 

The linearized 3D Euler equation around a 3D shear (U (y, z) , 0, 0) is 

(4.1) d t u + Uu x + vU y + wU z = -P x , 

(4.2) dtv + Uv x = -P y , d t w + Uw x = -P z , 

(4.3) u x + v y + w z = 0, 

with periodic boundary conditions. There are almost no results about the the linear insta- 
bility of general 3D shears. So we construct unstable 3D shears near unstable 2D shear flows 
{Uo (y) > 0, 0) where U (y) is periodic with period L y . First, we give a sufficient condition for 
linear instability of 2D periodic shears, which generalizes the result in |L03j for shear flows 
in a channel with rigid walls. 

Lemma 4.1. Consider a periodic shear profile U(y) G C 2 (0, L y ) with only one inflection 
value U s and 

(4 ' 4) ^ ^ -vwk > °- 

Let — a^ax ^ e ^ e l° wes t eigenvalue of the Sturm-Liouville operator 

(4.5) Lip = -ip" - K(y)ip 

with the periodic boundary conditions on y 6 [0, L y \. Then the Rayleigh equation 

(4.6) U"<f> -(U-c) U" - a 2 (j)) = 0, 



with periodic boundary conditions on y € [0, L y ] has unstable eigenmodes (Imc > 0) for any 
a G (0,a max ). 

Remark. Under the assumptions in the above Lemma, the lowest eigenvalue of L is negative, 
since (L (1) , 1) = — J K(y)dy < 0. A typical example satisfying ^4-4\) ^ s Uq (y) = sin ( j^y 



for which K (y) 



2 



Proof. The proof is similar to the case of rigid walls ( |L03j ). so we only point out some 
small modifications. Let cj) s be the eigenfunction of L corresponding to the lowest eigenvalue 
— a max- Then (0, c, a) = (cj) s ,U s ,a max ) is a neutral solution to the Rayleigh equation f)4.6p . 
By Sturm-Liouville theory, — a^ ax is a simple eigenvalue and we can take 4> s > 0. First, we 
study bifurcation of unstable modes near the neutral mode. Denote y\ to be a minimum 
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point of (p s and let y 2 — y% + L y . We normalize S such that <f> s (yi) = 1, <fi' s (yi) = 0. Define 
0i (y; c, e) and 2 ($/; c, e) to be the solutions of 

(4-7) - 0" + v ^ a _ c <t> + («Lx + e) = 0, 

with 0x (yi) = 1, 0' x = and 2 = 0, (f>' 2 (yi) = 1. Here e < and Imc > 0. Define 

J (c, e) = 0i (y 2 ; c, e) + 0' 2 (y 2 ; c, e) - 2 , 
then the existence of a solution to the Rayleigh equation (14. 7p with periodic boundary 
conditions on y G [1/1,2/2] is equivalent to the existence of a root of / with Imc > 0. When 
c — > 0, e — > 0— and |Rec| / Imc remains bounded, as in |L03j we can show that 0i (y;c, e) 
(02 (y;c,e)) converges to <f> a (y) (<f> z (y)) uniformly in C l [y u y 2 ). Here, <f> z (y) G C l [y u y 2 ) 
satisfies that 0' 2 (y 2 ) = 1 and Z (y 2 ) 7^ since Z (y) can not be another eigenfunction 
associated with the simple eigenvalue — « max . By similar calculations as in |L03| . it can be 
shown that when c — > 0, e — >• 0— and |Rec| / Imc remains bounded, 

-> 0^ (2/2) y 0^ (y) dy, 

and 

^ -0, (2/2) U g (|C/r ^) l,=a fe + V ^ (K (y) <Pl (y)) I (U (y) - U s ) dy 

Here, ai, • • • , a; are the inflection points such that U (a*-) = C/ s , = 1, • • • , Z and P J^ 2 
denotes the Cauchy principal part. Then by a variant of implicit function theorem as in 
[L03J . there exists sq < such that for any £0 < s < . there is an unstable solution e with 
c = c (e) to Rayleigh's equation (14.71) . By the same arguments in |L03j . such unstable modes 
can be continuated to all wave numbers a G (0, a max ). □ 

Our second step is to show that 3D shears near an unstable 2D shear are also linearly 
unstable. More precisely, we have 

Lemma 4.2. Let Uq (y) G C 2 (0, L y ) be such that the Rayleigh equation ( fj76| ) has an unstable 
solution with (ao,co) (ao,Imco > 0). Fixed L z > 0, consider U (y,z) G C 1 ((0,L y ) x (0,L 2 )) 
which is L y , L z -periodic iny and z respectively. If\\U (y, z) — Uq (y)\\ w i, P ^ l v )x(o l z )) (p > 2 ) ^ s 
small enough, then there exists an unstable solution e ta °( x ~ ct ^ (u,v,w, P) (y, z) to the lin- 
earized equation d^.i[ )-( f^Tg| ) with \c — cq\ small. Moreover, ifU (y, z) G C°° , then (u, v, w, P) G 
C°°. 

The proof of above lemma is almost the same as in the case of rigid walls ( [LLllj ). so 
we skip it here. By Lemmas 14.11 and 14.21 there exist linearly unstable 3D shears Uq = 
(U (y,z) ,0,0). Below, we show that the assumption (A2) of linear exponential dichotomy 
is true in spaces W^ er = H m (m > 1 is integer) for such unstable 3D shears. Then the 
assumption (A3) is automatic since fio = 0. Let Gt = e Lt be the linearized Euler semigroup 
near a steady flow uq (x) and denote r ess (Gt, H m ) to be the essential spectrum radius of G t 
in space H m . By rather standard semigroup theory (see e.g. [Shi83l Section 1 ]), to get the 
linear exponential dichotomy (A2), it suffices to show that r ess (G t ; H m ) = 1. The essential 
spectrum of linearized Euler operator had been studied a lot ( |FV91j |LM91j [SL09j |V96j ) 
by using the geometric optics method. We use the following characterization of r ess (G t ; H m ) 
in [5L09| . 
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Lemma 4.3. |SL09j Consider the following ODE system 
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(4.8) 



X t = Uq (x) 

£t = -du {x) T <f 
b t = —du Q (x) b + 2 ( du (x) b, £) £ 



where uq (x) is a steady flow of 3D Euler equation in T 3 and x G T 3 , £, b G R 3 . Denote 



(4.9) 



A m = lim - In sup 

6o-L|o, bo =1 



where [x (t) ,b(t) , £ (t)J zs i/ie solution of Hj4.S\) with initial data \Xq, £o, &oJ ■ 27ien we /lave 

r ess (G 4 ;if m ) = e' A "\ 
Lemma 4.4. For u = {U (y, z) , 0, 0) in T 3 , we have A m = 0. 
Proof. Denote 

x(t) = (x (t) ,y(t),z (t)) , i(t) = (6 (t) , £ 2 (t) , £ 3 (*)) , b(t) = (h (t) , b 2 (t) , b 3 (t)) , 

and 

% = (*o, l/o, *b) , fo = e 2 °, £ 3 °) > &o = (&?, 6°, 6°) . 
The solution of first two equations of (14.81) yield 

x(t) = xq + U (y , Zq) t, y (t) = y , z (t) = z 

and 

6 (t) = £?, 6 (*) = -tfytf* + £ 2 °, 6 (t) = -u g $t + £ 3 °. 

Plugging above forms into the equation of b (t), we have 

2 (£i°) 2 (Uyb 2 + U z b 3 ) 



(4.10) 



6l = - (Uyh + U z b 3 



at) 



(4.11) 



2^ (Uyb 2 + U Z b 3 ) - Uy^t) 



at) 



(4.12) 



2^ (Uyh + (^3° - Uzgt) 

at) 2 
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To show that A m = 0, it suffices to prove that 



b(t) 



only has polynomial growth, uniformly 



in (x ,£ ,b j . From equations (14. lip and (14.121) . we have 

J t (Uyb 2 + U z b Z ) 



2 (Uyb 2 + U Z b 3 ) 


'mUy + &$U M - (tf f (ui + m) t 




at) 


2 



(Uyb 2 + U Z b 3 ) d 



at) 



dt 



at) 



by noting that 



at)\ = fe ) 2 + {-Uy&t + e 2 ) 2 + (-u z & + 
= 1-2 ($&u v + mu z ) t + (e?) 2 (u 2 y + ui) e. 



Thus 

and 
(4.13) 



d 
dt 



(V y b 2 + V z b 3 ) £(t) =0 
U y b° 2 + U z bl 



(Uyb 2 + U Z b 3 ) (t) 



For any fixed t > 0, we find a lower bound for 



at) 



at) 



by minimizing the function 



/ (& e 2 °, e 3 ) = fe ) 2 + (-Uy&t + e 2 ) 2 + (-u^t + e 3 ) 2 

subject to the constraint (^?) 2 + (£ 2 ) 2 + (£!s) 2 = 1- By calculations of Lagrange multiplier, 

we get 



mm 

\i \=i 



at) 



2 2 + (U 2 y + Ul) - J (2 + (VI + vt) t*y - 4 



2 
2 



2 + (0* + 02) t 2 + ^( 2 + (pi + 0J) _ 4 
So for t > 0, we get the estimate 

1 



at) 



> 



2+(U* + U*)t*' 



Thus by KIM . 

\(V y b 2 + V z b 3 ) (t)\ < Cl t 2 + c 2 , 
for ci,c 2 > independent of (i?o,£o,^o)- By ( 14. 10|) -( l47T2|) . we have 



<2\{Vyb 2 + V Z b 3 ) 001, 



< \(Uyb 2 + V Z b 3 ) (t)\ 
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and thus 



b(t) 



< c 3 t s + c 4 



for some constants c 3 , c 4 > 0. This finishes the proof of the lemma. 



□ 
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